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Abstract 

We compute various correlation functions at the planar level in a simple supersymmetric matrix 
model, whose scalar potential is in shape of a double-well. The model has infinitely degener- 
ate vacua parametrized by filling fractions v± representing the numbers of matrix eigenvalues 
around the two minima of the double-well. The computation is done for general filling fractions 
corresponding to general two-cut solutions for the eigenvalue distribution. The model is mapped 
to the O(n) model on a random surface with n = — 2, and some sector of the model is described 
by two-dimensional quantum gravity with c = —2 matter or (2,1) minimal string theory. For the 
other sector in which such description is not possible, we find new critical behavior of powers of 
logarithm for correlation functions. We regard the matrix model as a supersymmetric analog of 
the Penner model, and discuss correspondence of the matrix model to two-dimensional type IIA 
superstring theory from the viewpoint of symmetry and spectrum. In particular, single-trace 
operators in the matrix model are naturally interpreted as vertex operators in the type IIA 
theory. Also, the result of the correlation functions implies that the corresponding type IIA 
theory has a nontrivial Ramond-Ramond background. 



1 Introduction 



Matrix models given by dimensional reduction of supersymmetric Yang-Mills theories have 
been attractive frameworks capturing nonperturbative aspects of superstring/M theory, 
which are based on conjectures [HO [3] and on AdS/CFT correspondence [HE]. Although 
such matrix models have been vigorously investigated since these proposals, it is still hard 
to unveil nonperturbative string vacua and understand their essential properties. 

In this situation, we consider a simple supersymmetric matrix model discussed in [6], 
which is a zero-dimensional matrix model analogous to two-dimensional Witten's model 
containing a real scalar field and a Majorana fermion [?]. Its scalar potential takes the 
form of a double-well. This model seems much more tractable compared to the matrix 
models of supersymmetric Yang-Mills type, and we hope the analysis here is helpful to get 
insights into nonperturbative dynamics of strings not ever clarified and in particular cru- 
cially related to target-space supersymmetry. Ref. [6] found that the model has infinitely 
degenerate supersymmetric vacua at the planar level. The degeneracy is parametrized by 
filling fraction v + (z/_) corresponding to the number of eigenvalues of a scalar matrix 
around the right minimum (the left minimum) of the double-well. The support of the 
eigenvalue distribution is two intervals for general v± (two-cut solutions). 

In this paper, we compute various correlation functions of the matrix model, and inves- 
tigate their critical behavior as the two intervals of the eigenvalue distribution approach 
to touch each other. Although this model is equivalent to a matrix model for the 0(n) 
model on a random surface with n = —2 [SJ |H1 QUI EH], analysis ever done in the literature 
is for one-cut solutions Q, where the eigenvalue distribution has the support of a single 
interval. Thus, the analysis here would be also helpful to consider a multi-cut case in 
general 0(n) matrix models. As a result of the analysis, we find critical behavior of pow- 
ers of logarithm in planar correlation functions of odd powers of the scalar matrix and of 
fermionic matrices. Such behavior has not ever been seen in the literature, although it is 
somewhat reminiscent of the logarithmic scaling violation in a matrix quantum mechanics 
called the c = 1 matrix model [13 EH El [IE] |. 

The Penner model [20] and its extension including source terms for "tachyons" (the 
Kontsevich-Penner model) [2TI 122] are zero-dimensional matrix models whose critical be- 
havior is the same as the c = 1 matrix model 0. They describe bosonic string theory 
in two-dimensional target space: (Liouville direction) x (S 1 with the self-dual radius). It 
suggests us an intriguing direction to regard our supersymmetric matrix model as a su- 
persymmetric version of the Penner model and to consider correspondence of the matrix 
model to two-dimensional string theory with target-space supersymmetry. Indeed, two- 
dimensional type II superstring theory with the same target space, (Liouville) x (S* 1 with 
the self-dual radius), is constructed [2U [251 ESI EZ] , where the target-space supersymmetry 
can exist only at the self-dual radius of the circle. 

1 For example, see refs. [T2"l IT31 RH] . 
2 For a review, see ref. [19] . 

3 Also is the normal matrix model [23], which corresponds to c = 1 noncritical strings on S 1 with a 
general radius. 
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The rest of this paper is organized as follows. After a short review of the supersym- 
metric double- well matrix model in the next section, we compute the one-point, two-point 
and three-point correlation functions among operators of the scalar matrix -^tr <p n at the 
planar level, and investigate their critical behavior in sections El H] and El respectively. 
The two-point functions of fermionic matrices are calculated in section |6j We summarize 
the results obtained so far, and discuss correspondence of the matrix model and two- 
dimensional type IIA superstring theory from the viewpoint of symmetry and spectrum 
in section [7J Analysis of Schwinger- Dyson equations for the matrix model is presented in 
appendix [A] We discuss dependence on the filling fractions of correlation functions among 
operators of the scalar matrix in appendix [Bl Some useful formulas for the calculation 
is given in appendix O 

The next paper [23] is devoted to computing various amplitudes in the type IIA theory, 
and to discuss the correspondence in detail by comparing the amplitudes in both sides of 
the matrix model and the type IIA theory. 

2 Supersymmetric double-well matrix model 

Let us start with a brief review of a supersymmetric matrix model discussed in ref. [6]. 
The action is 

\b 2 + iB(4> 2 - /i 2 ) + + 



S = Ntr 



(2.1) 



where B and <fi are N x N hermitian matrices, and ip and ip are N x N Grassmann-odd 
matrices. S is invariant under supersymmetry transformations generated by Q and Q: 

Q<p = ip, Qip = 0, QtP = -iB, QB = 0, (2.2) 

and 

Q<j> = -$, QV = 0, Qip = -iB, QB = 0, (2.3) 
from which they are nilpotent: Q 2 = Q 2 = 0. The partition function is defined by 



Z={-\Y J d, B d (f> d ip d ip c , (2.4) 
where we fix the normalization of the measure as 

d JV^ e -JVtr(|^) = J dN 2 Be _ NtT{ i B2) = ^ ( _ 1)jy2 J (^N 2 ^ d N 2 ^ e -Ntrm = L 

(2.5) 

We express the expectation value of a single trace operator -^tr O in the -^-expansion as 
l -*°) = (^0) +l^(i*o) +i( / ^ro\ +•--. (2.6) 



N I \N / N 2 \N / 1 iV 4 \N 
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(•) h is an ^-independent quantity representing an amplitude over random surfaces with 
handle h. Similarly, the expansion of the connected correlator of K single-trace operators 
is 

f\h i0k ) =£]v^^(n^ tr0fc ) (2 - 7) 

=1 la h=o \fc=i / c ,h 

with {-) Ch ^-independent . 

After integrating out B, we find that the scalar potential of <fi is in shape of a double- 
well: |(0 2 — /i 2 ) 2 . As shown in ref. [B], there are infinitely many large- A" saddle points pre- 
serving supersymmetries for p 2 > 2. The eigenvalue distribution p(x) = (-^tr<5(x — </>)) 
has the support Q = [—b, —a] U [a, b] with 



-2 + /i 2 , b=^/2 + p 2 , (2.8) 



and its explicit form is 



/\_{ v x \/ (x 2 — a 2 )(b 2 — x 2 ) (a < x < b) ,^ g\ 

\ — \x\ \/ (x 2 — a 2 )(b 2 — x 2 ) (—b < x < —a). 

The filling fraction v + (z/_) expresses the amount of the eigenvalues around the right (left) 
minimum of the double-well, satisfying v± > 0, i> + + v_ = 1. The cases of (v + , z/_) = (1, 0) 
and (0, 1) reduce the support of p(x) to a single interval of [a, b] and [—6, —a], respectively. 
The large- A" planar free energy F Q , which is the leading term of 

1 00 i 

F ^-Jf2^ Z = Y,J^i F ^ ( 2 - 10 ) 

h=0 

calculated from (12. 9p vanishes, and all of (-^tr£>™) [n = 1,2, •••) do so. It indicates 
that the solution preserves the supersymmetry. Note that the end points of the intervals 
(12. 8p do not depend on v±. Because it is not the case of bosonic double- well matrix 
models [SHIED], it is considered to be a characteristic feature given in supersymmetrizing 
the double- well matrix model. It is easy to see from (12. 8p that the supersymmetric two-cut 
solution (12.91) ceases to exist for /i 2 < 2. In fact, in ref. [31] it is found that when /i 2 < 2 
there is only a vacuum without supersymmetry and that we have a third-order phase 
transition between supersymmetric and supersymmetry breaking phase at fi 2 = 2. Hence 
/i 2 = 2 is a critical point in our model (12. ip and we will elaborate on critical behavior of 
correlation functions in the supersymmetric phase as fi 2 — > 2 + 0. 

The model are interpreted as the 0(n) model on a random lattice with n = —2, whose 
critical behavior is described by the c = — 2 topological gravity [H] (the k = 1 case of 
(2, 2k — 1) minimal string theory constructed by the one-matrix model [32l [33l EI]). The 
string susceptibility exponent is 7 = —1. The partition function (12. 4p after B, ip and ip 
integrated out becomes a Gaussian one-matrix model by the Nicolai mapping H = (p 2 —fi 2 , 
where the integration is over region where H + fi 2 is a positive definite hermitian matrix, 
not over all the hermitian matrices. Ref. [351 discusses that since the difference of the 
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integration region has only effects which are nonperturbative in 1/N, the model can be 
regarded as the standard Gaussian matrix model integrated over all the hermitian matrices 
at each order of genus expansion. 

The Nicolai mapping changes the operators -^tr (f> 2n (n = 1, 2, • • • ) to regular operators 
■j^tr (H + fi 2 ) n with respect to H. Hence the behavior of their correlators is expected 
to be described by the Gaussian one-matrix model (the c = —2 topological gravity) 
at least perturbatively in 1/N. Actually, we see that ^tr0 2n ^ is a regular function 

of /i 2 in section [3J and that the planar expectation value of the resolvent itr 

obeys the semi-circle law from ( 1A.16|) in appendix [A] However, the operators -^tr 
(n = 0, 1, 2, • ■ • ) are mapped to ±^tr (H + /i 2 ) n+1 / 2 that are singular at H = —fx 2 . They 
are not observables in the c = —2 topological gravity, while they are natural observables 
as well as -^tr0 2n in the original setting (12. 4p . In the following sections, we will see that 
correlation functions among operators □ 

ltr0 2n+1 , ^tr^+\ ^tr^ 2n+1 (n = 0, 1, 2, • • • ) (2.11) 

exhibit logarithmic singular behavior of powers of ln(/i 2 — 2) at the planar topology. 
Note that it is not straightforward to obtain correlation functions of such operators from 
Schwinger-Dyson equations as seen in appendix [S] 



3 One-point functions 

The expectation value of -^tr0 n (n — 1, 2, • • • ) in the large- iV limit is expressed in term 
of the hypergeometric function as 

/^-tr0"\ = / dxx n p(x) = {v + ±uJ)- [ dxx n+1 ^{x 2 -a 2 ){b 2 -x 2 ) 

\ N / Jn ' 7T J a 

= (,/+ ± (~ ! 3; , (3.1) 

where the plus (minus) branch of (u + ± vS) is taken for n even (odd). Note that when 
n is even, the hypergeometric function F (— f , 1,3; 4) is reduced to a polynomial of li 2 
giving nonsingular behavior with respect to li 2 : 

On the other hand, when n is odd, it exhibits logarithmic singular behavior for u = 

4 Note that tiip 2n = tr?/; 2 " = (n — 1,2, • ■ •) from the cyclicity of the trace and the Grassmann 
nature of ip and tp. 
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Z/j- — V- 



v> — v_ 



1/.L. — V. 



64 16 



oj + -u 2 \iauj + O(oj 2 



15ir 3tt 7r 

1024 128 16 2 I .,, 3 
1 oj H w 2 + -ur In w + 0(ur 

1057T 57T 7T 7T 

8192 2048 128 , 160 , 10 4l m , 4 

77— + — — u + oj 2 + —u^ + — w 4 lnw + 0{uj 4 

Slbii 21ti 7i Sir it 



In general, for k — 0, 1, 2, 
1 



sing. 



2 fe + 2 (2& + 1)!! fc+2 



7T 



(A; + 2)! 



w In oj + (less singular) 



(3.3) 



(3.4) 



where | S i ng . stands for removing entire functions (polynomials) of oj. 

The singular term of In a; becomes dominant after taking /i 2 -derivatives (k + 2) times, 
namely 



Qk+2 



1 

— tr 



k 2k+l \ 



(2/c + 1)!! 



C,0 



1/4- - 



7T 



lnw + O{oj°, ojlnoj) 



(3.5) 



4 Two-point functions of bosons 

In computing higher-point correlators ^n^=i J?^ T 4> ni J & t the vacuum with general filling 

fractions it is useful to reduce them to those at the vacuum with {y + , uJ) = (1, 0). 

It is discussed in appendix [Bj where 

ik \ I k \ 

( n ) =(*+- ^ ( n ^ tr ^ ) 

\*=1 / C,0 \ i=1 I C,0 

is shown up to K = 3. Here, the superscripts (z/ + ,z/_) and (1,0) are put to clarify the 
filling fractions of the vacua at which the expectation values are evaluated, and Jj counts 
the number of odd integers in {iii, • • • , tik}- 



4.1 Eigenvalue distribution with source 

In order to obtain higher-point correlators of jftr <fi p (p — 1,2, • • •), we introduce source 
terms Yl^Li iptr 4^ to the partition function: 




In the large- N limit, the eigenvalue distribution Pj(x) satisfies the saddle point equation 



dypj(y) P 



+ P- 



X 



x + y 



x 



fl 2 X 



P =i 



Pip x p-i 



(4.3) 



Let us consider the case of the filling fractions (1,0) % with the support of Pj(x) [cij,bj] 
(0 < a.,- < bj). We change variables as 



x 



A + Br] with 



A 



B 



(4.4) 



and put p(r]) = j-pj(y), to simplify (14. 3 p as 



2.y 



1 

B 



dr)p(r])P- 

i Z-v 



P =i 



(4.5) 



for £ G [—1,1], where p is normalized by dripijf) = 1. 

We act §_ x d£ \J\ — £ 2 P^z| 1° both sides of ( I4.5p . and apply the formula (10.8|) in 
appendix O Then we obtain 



p(0 



2-(A-p 2 )BC-B 2 (C 



00 BP 1 1 



(4.6) 



after using f lC.ip and ( 1C.2j) . The condition p(( = ±1) = determines A and B as 

2B 



A 



B 



+ ^ Jp |(X + B)S- lF (-| + 1 ,i, 1; 

i+ £fl(l- i ) B2( - 4+B)! " F (-l +2 -l 3i 



25 
A + 5 



1/2 



(4.7) 



(4.8) 



5 Curiously, the saddle point equation (|4.3p is sometimes not consistent for j p ^ with p odd, when 
the support of Pj(x) consists of two intervals f2j = [— bp — a'j] U [a,j,bj] (0 < aj < bj, < a'j < bj). If 
there is a point x € flj such that also — x G fij, eq. (|4.3I) for a; and that for —x contradict each other, 
because the l.h.s. of (|4.3p is odd for x — > —x, while the r.h.s. is not due to the source j p (p: odd). 
The l.h.s. tells that the eigenvalue at the position x feels repulsive forces from the other eigenvalues as 
well as from the mirror images of all the eigenvalues. Since it becomes an odd function of x, total force 
is not balanced unless the matrix model potential is even. Thus, the deformation by the source j p (p: 
odd) makes the general two-cut solution (|2.9p unstable, and would drastically change the support of the 
eigenvalue distribution. 

To avoid this difficulty, we first solve (|4.3|) for the filling fractions (1,0), and then obtain the amplitudes 
for general (v + , vS) by using (|4.1[) . 
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from which A and B are obtained iteratively with respect to {j p }. Up to the first order 
of {j P }, 



+ o(f 

4 



(4.9) 



+ 0(f), (4.10) 



CO s "I A 

p=l x 

p=l x 

where C(j 2 ) means a quantity of the quadratic order of {j p }- 
4.2 Two-point amplitudes (^tr (f) p jjtr <7 ) c - o 

Let us express the planar expectation value of O under the partition function with the 
source terms (14. 2 j) as ■ The cylinder amplitude at the vacuum with the filling 

fractions (1,0) is given as 



(1,0) 



C,0 



dj P \W /o 



c?) 



{3p}=0 



9j 



- /'dC (A + B()lp(() 
p J-i 



(4.11) 



Op>=o 



Plugging (Q)l . ( BSD and ^^IQj 1 int o fl^TT]) leads to 



(1,0) 



C,0 



rfC <9 C [(/i 2 + 2C) « v 7 ! 3 ^ (<% A + d jp B C) 



i 



{j P }=o 



P 



2tt 



2 / ^c^4=^ [ 1 dav^ms+2^ 

J-i vi-C 

1 



xP- 



(4.12) 



Note that the linear-order contributions to A and B with respect to appear in the 
first term, but it is a total derivative term and vanishes. Thus, from ( 14. ip the amplitude 
at general filling fractions is 



1 




(f+>f-) 








iv tr 




C,0 


("+ 


; 2tt 2 


/> 



(/i 2 + 2C) 



f 



l-C 2 J -i 



^v / i 3 e(^ + 20 f_1 p 



(4.13) 



where (1 = 2 for p and q odd, (1 = for p and q even, and (1 = 1 for one of p and q odd. By 
integrating by parts (I4.13P with respect to £, we obtain a symmetric expression for p and 
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i i \ (y+,v-) 



W+. — V' 



2tt 2 



C,0 

'^(^ + 20 
1 



^^fp7^(a-D.(4.i4) 



i-c 2 V-i ^A 3 e (c-o s 



4.2.1 Case of p and g odd 

Let us compute further (14. 13|) for case of p and q odd. We can put p — 2k + 1, q — 2£ + 1 
with k < £, i.e. £ = k + m (m = 0, 1, 2, • ■ ■ ), without loss of generality. Symmetrizing 
(I4.13P with respect to ( and £ leads to 

Itr^+'ltr^+A 



/ 



C,0 



X 



2fc + 1 r 1 (/i 2 + 2C) 



47r 2 



fe-1/2 /•! 



l - c 2 



(/i 2 + 2£) 



fc-1/2 



-1 



i - e c - £ 



[Gu 2 + 2C) m+1 (l - £ 2 ) - + 2£) m+1 (i - C 2 )] • 



(4.15) 



Here and in what follows, we omit the superscript (z/ + ,z/_) when it does not cause any 
confusion. We express the amplitude in terms of 



v- 



n(A(l+uj)y-iF( i-pi, 1; ' 



1 + U) 



(p = 0,l,2, 



(4.16) 



as 



^C,0 
m 

.p=l 



i + (4 - 



m—p+k- 



(4.17) 



Note that I p _i yields logarithmic singular behavior (~ a; p lna;) for u; ~ +0. For p 7^ 0, 



however i_ 1 has a finite limit as u; — > +0. In fact we find 

f 2 



(On - IV 

_ 2P -i I ^ )■■ UJ P\ nu + O(co p+l In w, co p ) 
p\ 



+ ( 2 3p " 



(4.18) 



with (-1)!! = 1. 
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The explicit form of the first several amplitudes are 



^tr0itr0 



^tr0ltr0 3 



) 



/ c*,o 



C,0 



\ 



C,0 



N^N**/ 



\ 



tr J — tr^ 5 
N N 



sing. 



sing. 



sing. 



sing. 



2tt 2 



w(lnu;) 2 + C(w lnw) 



4 3 

oo In w - — w 2 (ln w) 2 + C(w 2 In oo) 



7i 

24 



2tt 2 
6 



— oo 2 In w -u; 3 (m w) 2 + C(w 3 In w) 



71 + 



C,0 



smg. 



\ 



C,0 



= (is + - V. 



smg. 



— - oo In oo H — -cu Inw cu Inw 

+ O{oo 3 lnw) 

32 144 o, 45 4/1 

— a; muH —a; In a; -oo (Incur 

7T 2 7T 2 27T 2 



+ 0(oo A lnw) 
640 



3, 720 4, 90 5/, 

cu In uH —oo in oo oo (ma; J 



3tt 2 
+ C(w 5 In a; 



71 + 



7T + 



(4.19) 



A characteristic feature of the amplitudes is singular behavior of (lncu) 2 multiplied 
by powers of oo. The leading singular term of (lncu) 2 behavior in (-^tr0 2A:+1 -^tr 2m ) co 
with £ = k + m (m — 0, 1, 2, • • • ) takes the form: 



C*,0 



=-(" + -"-) 2H ^ 2 * + ' 



(In w) 2 -leading 

(2p + 2A; - 1)!! (2m - 2p + 2 A; - 1)!! 



E 

.p= 
+2 



^ (p + A;)! (m-p + fc + l) 
(2A: — 1)!! (2m + 2ife- 1)!! 



fc! 



(m + k)\ 



oo 



2fc+m+l 



(lnw) 2 . 



(4.20) 
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4.2.2 Case of p odd and q even 

Let us compute the amplitude for case of p odd (p = 2k + l) and q even (q = 2£) in (14. 13ft : 



— ti(f) 2k+1 — tr^ 
N 



[V< - V_ 



2k + 1 f 1 (^ + 2(Y f 1 



Cfi 



2tt 5 



Use of 



and (1C.9I) gives 



^ + 20 £ = £ (f) (/i 2 + 2e)^ y(c - o r 

r=0 ^ ' 



(4.21) 



(4.22) 



■ i 

i v 7 !^ C-{ 



-EE r (P - 



;-2) 



r-2m 



r (£ — r)!(r — 1 — 2m)! (m!) 5 



(^ + 2o £ - r e 



r tr-l-2m 



(4.23) 



r=l m=0 

with [x] the greatest integer not exceeding x. Plugging this into (I4.2ip . we have 



C,0 



[ 2 ] r—l—2m 



:-2) 



r-2m+s+l 



1 r— . « C 1 

<"+ - ] '-}-( 2k + 1 )z2 J2 r (£ — r)\s\(r — I — 2m — s)\ (,„!)-' 

r=l m=0 s=0 V 7 V / V / 



x(4(l + w)) 



fc+£-r-i 



55 + - F -fc-Hr + -,s + -,s + 3; 



1 



1 + 
(4.24) 



from which logarithmic singular behavior of the amplitude is seen as 



tr0 2fe+1 — ticf) 



N 



C,0 



sing. 



2^ (2fc + l)!!(2l-3)!! k+l 
2tt (Jfe + 1)! (£-l)\ 

+ (less singular) . (4-25) 
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The first several amplitudes are explicitly given as 



iv^ tr ' 



^tr0ltr0 4 



1 n 1 j2 

iv tr0 V r0 



^tr0 3 ltr0 4 



1 ,5 1 ,2 



^tr0 5 ltr0 4 



C,0 



C,0 



C,0 



C,0 



C,0 



sing. 



sing. 



sing. 



sing. 



sing. 



sing. 



V- 



1 3 

— w In co + — co 2 In co + O(co 3 In a;) 

7T 07T 



— In co H cj 2 In co + (9(u; 3 In co) 

7T 47T ^ 



3 3 

- to 2 In w + — to 3 In w + O(co 4 In a;) 

7T 47T 7 



6 39 

- co 2 In co + — a; 3 In co + O(co 4 In w) 

7T Z7T 



— cu 3 lna; + — a/lnw + 0(u 5 lnu) 



— to 3 In co + — cu 4 In w + O(co 5 In 

7T 47T 



(4.26) 



4.2.3 Case of p and g even 



The amplitude (14.131) becomes a polynomial of co for case of p and g even. In this case 
(I4.13P can be derived independently by (IA.21I) from which we read off 



N 



:tr< 



N 



:tr« 



1, 



C,0 



AT 



:tr« 



iV 



:tr. 



4(1 + 2w), 



C,0 



(4.27) 



4.3 Operator mixing 

In this subsection, we discuss operator mixing which makes the logarithmic singular be- 
havior of the cylinder amplitudes more transparent. Eq. (14. 25p indicates that the leading 
logarithmic behavior co k+1 In co is brought by the operator of an odd power -^tr 2fc+1 since 
it does not have ^-dependence. (The operator of an even power -^tr <p 2i just affects the co- 
efficient of the u k+l Inu term.) On the other hand, from (I4.19P we find that in the cylinder 
amplitude of operators of odd powers of 0, the (In a;) 2 terms of our interest are in general 
less singular than some In co ones. These observation suggests that we can construct a set 
of new operators by mixing each operator of an odd power and operators of lower even 
powers in such a way that the leading singular behavior in the cylinder amplitudes of 
the new operators will be given by (lncu) 2 terms. Note that operators of even powers we 
add do not affect (In a;) 2 terms at all, because cylinder amplitudes involving even-power 
operators do not have a singularity of (In a;) 2 , as seen from (I4.25P and section 14.2.31 
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4.3.1 Explicit construction 



It is easy to check that this prescription works at least lower orders by explicit construc- 
tion. Since the most singular term in (^tr <p j^ti 0) co is given by co(lnuj) 2 as in (I4.19p . 
we do not have to do anything. Let us set 

$i = ^tr0. (4.28) 

Then (14.191) tells us that ($i-^tr0 3 ) has the ulnu term which is larger than u 2 {lnuj) 2 . 
Thus we combine 4ftr0 3 with an operator of the power which is even and smaller than 



as 



$3 = ^tr< 



U>4 



z/_)c^tr 2 , 



(4.29) 



where c is a constant fixed by imposing that singular terms of ($1^3) start from u 2 (lnu) 2 . 
From (I4.19P and (14. 26p . c can be determined as c = — - so that 



($i$ 3 ) 



Wx. — V- 



sing. 



2tt 



lu 2 (Into) 2 + (less singular) 



(4.30) 



It is interesting that once we fix c = — -, it also absorbs the unwanted a; 2 In a; term in 
<itr0 3 ^tr0 3 ) CQ in (EH} as 



($3$3> Ci0 



sing. 



6 

— tu 3 (\nuj) 2 + (less singular) 

TV 2 



Likewise, we assume 



1 



1 



tr </> 5 + - a— tr 4 + - uJ) /3—tr <j) 



(4.31) 



(4.32) 



and fix a and /3 so that the leading singular term in ($£$5)^ (A; = 1,3,5) will contain 
(lnw) 2 . It turns out to be possible, if we allow w-dependence of a and j3 as polynomials 
of co of lower degree (entire functions of uS) as 



$1 = — tr ( 

1 N 



$ 3 
$5 



^tr0 3 + - (ag + ago; + 0(c 2 )) ^tr0 2 



iV 
1 



U>4 



5.2 '5,2 ^ • 0(w )J y" ' 



(4.33) 
(4.34) 

(4.35) 



with o4°2 = — ~. Notice that the even-power operators do not affect the singular behavior 

1 



of the one-point function in (13. 4ft : 

(*2fc+l) 



sing. 



N 
12 



:tr 



(4.36) 



sing. 



and that inclusion of a^\u + 0(u 2 ) in (14.34p does not spoil (I4.30p and (I4.3ip because 
it only affects higher order, hence less singular terms. Then vanishing of a; In a; term in 
($ 1 $ 5 ) C0 , of u 2 \nuj term in ($3$ 5 ) c , , and of w 3 lnc<; term in (& 5 & 5 ) C0 lead to exactly 
the same equation: 

= ^ + 2^ + A (4.37) 



Furthermore requirement for u 2 \yiuj term in (<&i$5) C0 , w 3 lnw term in ($3$5) c<0 and 
a; 4 lna; term in {§ 5 § 5 ) CQ to vanish give two independent equations (it turns out that all 
of the three are not independent): 



^ + HiL<°) + £ a (°) + 2 « (1) 1 <-' J ' 

77- 4 5,4 8 ' 2 5 ' 4 
104 39 (o) 3 (m a , , .n 



0=- + ^ + ^S + 24:i + 4i (4-38) 



= — + -agi + -a$ + 6a$ + 3a$. (4.39) 



From these equations, we see that the linear combinations 
$1 = — tr0, 



$3 = ^tr0 3 



,_)i(l + ogc + 0(c 2 )) Itr^ 2 . 
$ 5 = ltr0 5 - (u + - vJ) 1 (l + a$ u, + 0( W 2 )) ^tr0 4 

- (*+ " ^ (l + 3(1 - agj) w + O(co 2 )) itr 2 (4.40) 

with 0:33 = a^/tts^) ^54 = ^l/^l give the behavior 

($ 2fc+1 $ 2 m)co . = (f+ - V-Y [u k ,ico k+e+1 (\nu) 2 + (less singular)] (4.41) 



sing. 



for k, £ — 0, 1,2. Mjt ^ is a constant equal to the coefficient of the u k+e+1 (lnu) 2 term in 
(■^tr0 2fc+1 j^tr (f) 2£+1 ) c Q . an d ^54 remain undetermined from the requirements with 
respect to k, I — 0, 1, 2. (They could be fixed by considering amplitudes for higher k and 



4.3.2 General structure of operator mixing 

As we have seen above, the equations for afl, ■ ■ ■ , from the requirement in ($1 $ 5 ) co , 
($ 3 $ 5 ) C0 and ($ 5 $5) c are highly degenerate. Let us discuss why we have such degener- 
acy in general case. Recall that our prescription is to determine the polynomials a 2 £ + i j2 i(w) 
in 

1 1 1 

$ 2m = ^tr0 2£+1 + (u + - v_) ^« 2m 2i ( w )-tr0 21 , (4.42) 

i=i 
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so that f!4.4ip will be satisfied for < k < £. Since (-^tr 2fc+1 -^tr <P 2£+1 ) C0 contained 
in ( < f ) 2fc+i < ^ ) 2£+i)(7o P oss ibly have u k+1 lnu, u h+e lnu terms which are larger than 
uj k+e+1 (In uj) 2 of our interest, we have to tune the polynomials c^+i^i^) i n such a way 
that these £ terms all vanish for any < k < t. We show that for each fixed £, vanishing 
of the leading \nuj term (uj k+1 In uj) gives exactly the same equation for any k. Under 
(fl~4"2|) we have 



($ 2fe+ i$ 2 m)c,o = ( ^ tr ^ +1 ^ tr ^ +1 ) co 



i=l 



ll 

N' 



tr 



; 2V 



C,0 



+ (v + - uj) ^a 2 fc+i,2j(w) ^tr0 2j -^tr0 2£+1 ^ + (regular terms), 

(4.43) 



where the regular terms come from cylinder amplitudes of even-power operators. Let us 
concentrate on the leading In a; contribution in this equation. From (I4.17P and ( I4.18p . the 
leading In a; contribution reads for k < £ 



l tT(j) 2k + l]_ t ,2i + l 



N 



N 



In uj- leading 



2^ (,_!), (2^ + 1)!! 
' 2tt 2 (2£-l)\\ (Jfe + l)! ' 

(4.44) 



while for k 
1 



_ tr ^2fc+l^ t ,2/c+l 



N 



N 



o4fc or, I i 

- 2 ' ™ fc+1 lnu;. (4.45) 



In tu-leading 



7T 2 fc(fc+ 1) 



Thus from (14.251) the first and second terms on the r.h.s. in (14. 43p have the leading In a; 
term as uj +1 In uj, while the third term provides uj e+1 In uj. Hence our requirement amounts 
to for k < £ 



() = / i_ tr ^2*+ll r 



N 



N 



1=1 



C,0 
1 



In cu-leading 



N 



;tr< 



C,0 



(4.46) 



In w-leading 
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and for k = £ 



In w-leading 

(0) / 1 .10 i i 1 
i=i J ' / O,0 



+ 2(1/+ - !/_) iV^^V 1 



(4.47) 

In w-leading 



Important observation here is that the ratio between two cylinder amplitudes appearing 
these equations is independent of k. Namely, from (14.251) and (I4.44p we have for < k < £, 
\<%<i 



£tr0 2 ^tr^+% o 



from which (I4.46P becomes 



ln^leading^ K _^g2(^-l)! (i-l)! (44g) 



7T 



In w-leading 



(2£-l)!!(2i-3)!! 



i=l 



(i-1)! 2m - 2 ' 7T (2<-l)!!' 11 



The same equation is obtained for the case k — £. In fact, (I4.49P with £ = 1 reproduces 
-- obtained from (I4.30p and (I4.3ip . and for £ = 2 it gives the same equation as 



(o) 

"3 2 — 

(Q7|l . 

To summarize, we have identified the reason why we get the same equation for several 
k for fixed £ as in (14.371) . It originates from the universal structure of the leading In a; 
term that the ratio (I4.48P does not depend on k. We expect such a universal structure 
persists even in higher-order In a; terms as seen in (I4.38P and (I4.39p . These observations 
suggest us that it would be possible to determine a^+i^i^) (1 — * — ^) consistently 
by our requirement (I4.4ip . We emphasize that dependence on k and that on £ are not 
factorized in the (In a;) 2 term (14.201) . which is clearly different from the situation in the 
leading In a; term (I4.44j) . 

The prescription explained in this subsection is reminiscent of the universal and 
nonuniversal parts of amplitudes in matrix models for the two-dimensional quantum 
gravity. There, disk and cylinder amplitudes have regular terms with respect to the 
cosmological constant t, which are identified as contribution from surfaces constructed by 
small number of triangles in the dynamical triangulation. Such terms are nonuniversal, 
but larger than universal terms which are relevant to the continuum physics. Thus, in 
order to take the continuum limit, we have to subtract the nonuniversal parts from the 
amplitudes in advance [36J. 
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5 Three-point functions of bosons 

Similar procedure to the case of the two-point functions can be used in computing three- 
point correlation functions. For the vacuum with the filling fractions (1,0), 



itr^ltr^itr, 



(1.0) 
CO 



d 2 



djpdj 



p u Jq J -I 



daA + BQ^p{Q 



d 2 I 1 \ {j) 



{j P }=o 



{j P }=o 



(5.1) 



It turns out that the second-order derivatives of A and B with respect to {j p } appear in 
total derivative terms and vanish. In terms of J's in (I4.16p . we have 



Itr^itr^itr, 



(1.0) 



Co 



— r 
2t7 



d jp Ad jq B h__ x + -(d Jp Ad Jq A + d Jp Bd 3q B) (Ir - ^j fi _ x ) 



+ 



p + r 



d lq B+[d ]q A-^d ]q B^ 







{j P }=o 



{j P }=0 

1 1 \ (1 ' 0) 



+(p v), 



(5.2) 



y 3p J,x \{j }=o anc ^ ®ip^\{j }=o can ^ e rea d fr° m (14. 9 P and (14.101) . 



where djA\ 

The amplitude for general filling fractions (^tr (jf -jj^tr (j) q jjtr<f) r )q^' U ~^ is obtained by 
multiplying the expression (15. 2 j) by the factor (z/ + — as in (14. ip . 



Case of p = q 



1 . From 



d h A \{j p }=0 



d h B \{ jp }=0 



Flii.l; 



1 



' In a; + - ln2 + O(u\nco) 



4tt 



71 



. ,'3 3 
F -,-,3; 



1 



32(1 + w) 3 / 2 V 2 ' 2 ' ' 1 + UJ 
1 1 
= — lnw + -(l-ln2) + 0(wlna;), 

47T 7T 

the amplitude for the simplest case p = g = r = lis seen to behave as 

3^ 



1 

N' 



:tr< 



W+ — V- 



C0 



16tt 3 



(lnw) 3 + C((lno;) 



(5.3) 



(5.4) 



(5.5) 



when bj ~ +0. 
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Case of p = q = 1 and r = 3. Similarly, we see that the amplitude for p = q = 1 and 
r = 3 has the behavior 



N 



-tr< 



1 

TV' 



-tr< 



•— (lnu;) 2 + — (21n2- 1) lnu; 

2"/T 7T d 



+ - 



7T" 



(In2) 2 + hi2 + - + — u(\nu) 3 + 0(w(lnw) 



. (5.6) 



The terms in the first line in the r.h.s. look nonuniversal singular ones, and are larger 
than the term ^3 cu(lna;) 3 which seems universal. As is seen shortly, they are absorbed 
by the operator mixing in section 14.31 



Case ofp — q — l and r = 2. We obtain the behavior of the amplitude for p = q = 1 
and r = 2: 



1 



N 



1 



— trra — tr 1 



N 



C,0 



1 • ^ 2 I 



<,7T Z 



(\nu) 2 + — (21n2- 1) lnu; 

27T 2 



+ — (In 2 - (ln2) 2 ) + 0(u(\nuf 



7T 



.(5.7) 



Operator mixing. From (15.61) and (15. 7p . we see that the operator mixing (14.401) dis- 
cussed in section 14.31 absorbs the nonuniversal singular terms and leads to the desired 
result 

2 +^-io(\mo) 3 + O(oo(\noo) 2 ) . (5i 



7T" 



Higher-point functions. The results obtained for the one-, two- and three-point func- 
tions of operators <3>2fc+i (k = 0, 1, 2, • • • ) naturally suggest the form of higher-point func- 
tions as 



n 

ii=l 



2fc;+l 



C,0 



- z/_) n (const.) ^-T+EIUto-i) (i na; )n + ( less singular) (5.9) 



with 7 = —1. Besides the power of logarithm (\nu) n , it has the standard scaling behavior 
with the string susceptibility 7 = — 1 (the same as in the c = —2 topological gravity) and 
the gravitational scaling dimension k of $2fc+i ; if we identify u with "the cosmological 
constant" coupled to the lowest dimensional operator on a random surface [33, |38l [39] . 
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6 Two-point functions of fermions 

In this section, we consider the two-point function of fermionic operators 

_L tr ^ fe+ i_L tr ^2 fe+ i\ {k = o, 1, 2, • • • )• (6.1) 

It is convenient to diagonalize the matrix as = C/diag(Ai, • • • , A n (U G 
SU(N)) and to write 

B = UB'U\ ip = Uip'U\ if) = Uip'Ul (6.2) 

Then, 

Eq. (16.11) can be computed by using ( I6.3P as Wick contraction. The leading contribution 
in the large- N limit comes from planar diagrams: 



( 2k + l )Tfkr £ 



jnp2k-\ 

U,*2,— i*2fc+l 



i \ Kl + ^2 ^2 + ^3 ^2fc+l + / 

= (2k + l) / dxidx 2 ■ ■ • dx 2k+ \ p(xi) p(x 2 ) ■ ■ ■ p(x 2k +i) 
Jn 

xP P P . (6.4) 

Xi + X 2 X 2 + X 3 X 2k +1 + X\ 

6.1 (^tr^tr^ 

The k = case in ( 16. 4p is given by a disk amplitude (13. ip with n = — 1: 
-tr^-tr^ = - 2 J n dx-p{x) 



2 v v 77 V 2 2 1 + w 

4 1 

h -wlncu + O(w) 

3n ix 



u ~ +0), (6.5) 



exhibiting In singular behavior. Supersymmetry invariance implies that this is equal to 
(-^tr(iB) -^tr0) co = \-§^ (^tr0) Q , interestingly which can be seen from ( 13. 3p . 
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6.2 (itr^ 3 ^tr^) C0 

We should treat carefully the product of the principal values in (16. 4p for k > 1. Let us 
compute it in the k = 1 case: 

^ M 4 M7 } G0 = 3 / fi ^ dz "( I »"«' , W p ^ p ?TI p J^' (6 ' 6) 

We write the principal values as 

1 1 x — y + i5\E 



P 



x + y 2 2 -— ' x 2 — y 2 — e 2 + i5ie2x' 

8i=± 

1 1 y — z + z5 2 e 

y + z 2 ^— ' y 2 — z 2 — e 2 + iS 2 e2y 1 

02=± 

P^- = ly + (6 . 7) 

z + x 2 z 2 - x 2 - e 2 + iS 3 e2z 1 ; 

with e(> 0) sufficiently small. Using 

(x - y + i5 1 e)(y - z + iS 2 e)(z - x + iS 3 e) 

= z(x 2 — y 2 ) + i5ie(y — z)(z — x) + (cyclic permutations) + 0(e 2 ) 
= z(x 2 — y 2 — e 2 + i5\e2x) + i8\e[— (y + z)x + (y — z)z] 
+ (cyclic permutations) + 0(e 2 ), (6.8) 

where "cyclic permutations" means the simultaneous cyclic permutations of (x, y, z) and 
(Si, 62, S3), we obtain 



x + y y + z z + x 

1 1 



1 y 

2 3 

<5l,<52,53=± 



y 2 — z 2 — e 2 + iS 2 e2y z 2 — x 2 — e 2 + i<53e2,s 
i<$ie{-(y + z)x + (y - 2)2} 1 



x 2 — y 2 — e 2 + %Sie2x y 2 — z 2 — e 2 + iS 2 e2y 
1 



z 2 — x 2 — e 2 + i(5se2z 
+ (cyclic permutations) + 0(e 2 ) , (6.9) 
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and thus from (16. 6p 



1 „1 73 

— trw a — tn/r 
N N 



3-3 
"23" 



^2 / dxdydz p(x)p(y)p(z) z 



+ 



y 2 — z 2 — e 2 + i5 2 e2y z 2 — x 2 — e 2 + i5 3 e2z 

^2 I dxdydz p(x)p(y)p(z){-(y + z)x + (y - z)z] 
?_ x l J n 



(6.10) 



2 3 

<Sl,<52,<53=± 

z5ie 1 1 

x 



x 2 — y 2 — e 2 + i5ie2x ?/ 2 — z 2 — e 2 + i5 2 e2y z 2 — x 2 — e 2 + i53e2z 

(6.11) 

As e — >■ 0, the first term (16.10P is calculated as 

« = 9jjzp(z)z^Jjxp(x)P^^^Jjyp(y)P-^^ 

dzp(z) z(z 2 - p 2 ) 2 , (6.12) 

by using the saddle point equation 

/ dzp(z)P - 1 = i(x 2 - /i 2 ) (x G fi). (6.13) 
x — z 2 

In the second term (16.111) . since 

5i=± 

the vicinity of x = ±y has major contribution to the x-integral. After the x-integral, the 
integrand will be proportional to e (p y2 l _ z 2 j ■ Then the z- integral could have nonvan- 
ishing contribution from the neighborhood of z = ±y, where the square of the principal 
value becomes singular. 

Let us evaluate (16.111) precisely by focusing the vicinity of x — ±y and z = ±y as 

x = ±y + ex, z = ±y + ez. (6.15) 

At x = y + ex and z = y + ez, the function 

A = y i6ie \ I ( 6 .i 6 ) 

. f-f , x 2 — y 2 — e 2 + i5±e2x y 2 — z 2 — e 2 + i5 2 e2y z 2 — x 2 — e 2 + i5 3 e2z 

Ol,O 2 ,03=± 
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becomes 

1—11 z z — X 

= ^~yT ^TfJ ~T^7 (z-x) 2 + l ( " ' 

We also get similar results for other cases in (I6.15p . 

Because the factor p(x)p(y)p(z){ — (y+z)x+(y— z)z} in the integrand varies sufficiently 
slowly compared with A around x = ±y and z = ±y, we may fix it to its value at x — ±y 
and z = ±y in evaluating the integral. Then, we obtain 

(EII]) = ^- / dyp(y)-(p(y) 2 -p(-y) 2 + 2p(y)p(-y)), (6.18) 

4 j(i y 



after the formula 



°° , , 1 ~Z ~Z-X TV 2 . 

dx dz————- — - — = — (6.19) 



x ,_ x ? + l? + l(i-i) 2 + l 3 
derived by picking up residues is used. Noting 

p{x)p{-x) = ^* 2 [4 - (x 2 - p 2 ) 2 ] = X ~ [V l~ V - ? A± ~ ^ ~ » 2 ) 2 } 
p{yf - p{-y) 2 = sg n(y)^±^y 2 [A - (y 2 - p 2 ) 2 }, 



(u + -u-) / dyp(y)\y\[A-(y 2 -p 2 ) 2 }= / dy p(y) y[A - (y 2 - p 2 ) 2 ] (6.20) 
Jn Jo. 

makes (I6.18P expressed in terms of a single p(y) as 

(EH]) = l[ dyp{y)y[A-{y 2 -p 2 ) 2 } 



-{u + -u-Y / dyp(y)y[A-(y 2 -fi 2 y\. (6.21) 



Thus, the final result is 



1 W^tn? ) = (6.12) + (6.21 ) 



N N 



Vitr [40 - 30(0 2 - p 2 ) 2 ] 



o 



- ^) 2 (itr [40 - 0(0 2 - ^ 2 ) 2 ] ^ . (6.22) 
It is a linear combination of the planar expectation values of odd-power operators com- 
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puted in section^ For uj ~ +0, (I6.22p behaves as 



1 ,3 1 73 



[Va. — V-\ 



O,0 



192 48 18 2l 9 , , „, 2 

+ uj u 2 In u (49 - 40 In 2) u 2 

35n 357r 7r 5ir 



117 



u 3 \nui + 0(u 2 



512 192 24 2 6 m/ 

1U57T 007T 07T 7T 



(6.23) 



We can check that this result can be obtained by a different treatment of the principle 
value from (16. 7p as 



P-LJy ! . 

x + y 2 f-^ x + y + iSie ' 

0i=± 



P J_ = Iy ! . 

y + z 2 f-^ y + z + iS 2 e ' 

02=± 



1 1 Z — X 

-1-1" 9 ^— ' — ^2 — 



2 — x + i5 3 e 



z + x 2 ' z 2 — x 2 — e 2 + z<5 3 e2,2 

<5 3 =± 



(6.24) 



This strongly supports the validity of the computations here. 



6.3 Operator mixing 

In the derivation of (16.231) . the factor (u + — uJ) 3 appears in a nontrivial way from the 
product of three p's in (I6.18p . For general k, it is expected that terms proportional to 
{y+ — V-) 2k+1 arise in (^tr ?p 2k+1 -^ti ip 2k+1 ') c Q from the product of (2k + 1) p's. Hence 

terms proportional to (v + — z/_) 2fc+1 are considered to be characteristic of the amplitude 
<itr^^tr^ fc+1 ) co . 

In (I6.23p . the leading singular term for u ~ +0 is not - (u + — u^) 3 u 3 \nu, but — — (u + — 
v_)u 2 1iiuj. Similarly to what was discussed for the operators $2fc+i m section |4TB| we can 
take a new basis of operators to make the (v + — vJ) 3 ui 3 In a; term dominant in the singular 
terms. Let us take the new basis as 

tfi = £tr^, *i = — tr$, 

^ 3 = itr tp 3 + (mixing), ^ 3 = — tr tp 3 + (mixing), (6.25) 

where "mixing" means operators to be added so that 

(*2fc+i*2<+i) ro =5 k/ v k (v + -v_) 2k+1 uj 2k+1 In lo + (less singular) (6.26) 

' sing. 

with Vf. constants holds for k, £ = 0, 1. 

As "mixing" operators in ^3, we consider fermionic operators which have lower powers 
of ip than -ip 3 and preserve the Q-exactness: ^3 = QE as ^tr-0 3 = Q (^tr (pip 2 ) . They are 
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allowed to be multiplied by polynomials of u as « 2 £+i,2i(w) m (14.421) . We find -^tr{(iB 
<p 2 + p? as an interesting candidate by noting that 



itr {{iB - 2 + /x 2 )^} itr {(iB - 2 + /i 2 )^} 



c.o 



f 

N' 



tr (i/ji/j) ) —— I dxdy p(x) p(y) P 



1/4. - V- 



1 

128 32 



— / dx p{x)x[x 2 — /i 2 ) 



4 



+ w + -w 2 lnw + — (23 - 24 In 2) w z + -w J lnw 



105tt 15tt 



7T 



3^ 



7T 



+—(67- 120 In2)w 3 + C(w 4 lnw) 

07T 



-o), 



(6.27) 



where (I6.3p and the saddle point equation were utilized in the second line, and that 



i_ tr{( , jB _ 2 + /i 2 m ^ tr ^ + l\ 



ItrV 2A+1 ^{(z^-^ + zi 2 )^} 



with = 0, 1. 

It can be shown from these equations that either of 



= 

0,0 
(6.28) 



and 



*3 

^3 

^3 
^3 



^tr ^ + JL (1 + u + (9( w 2 ))ltr {(.5 - 2 + 

Itr ^ 3 + JL (1 + w + 0(c 2 ))ltr {(iB - 2 + /i 2 )^} 



^tr ^ - ± (1 + co + 0(c 2 ))^tr {(zB - 2 + 
Itr ^ - A (1 + w + 0(c 2 ))ltr {(zB - 2 + 



(6.29) 



(6.30) 



does the job (I6.26P with v = - and v\ = -, where we have implicitly assumed the same 
coefficient in the mixing for ^3 and ^3. 

Eq. (16.261) indicates that \&2fc+i and ^fc+i have the gravitational scaling dimension k 
(the same as $2fc+i), besides the In a; correction. 



7 Correspondence to type IIA superstring theory 

We computed various amplitudes of the supersymmetric double- well matrix model (12. ip . 
and observed critical behavior with powers of logarithm in planar correlation functions 
among operators (12. lip , which are not belonging to the observables in the c = —2 topo- 
logical gravity. The result reminds us of the logarithmic scaling violation of the two- 
dimensional string (two-dimensional quantum gravity coupled to the c = 1 matter). Note 
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that the Penner model is known [20] as a zero- dimensional matrix model to describe the 
two-dimensional bosonic string whose target space is R x S 1 , where R is the Liouville 
direction, and S 1 is the matter direction compactified to the self-dual radius. Thus it is 
interesting to proceed with an interpretation of our matrix model as a supersymmetric 
version of the Penner model describing two-dimensional string theory with target-space 
supersymmetry. 

Let us begin with symmetry consideration. Suppose ip and ip are regarded as target- 
space fermions in the corresponding superstring theory, i.e. ip is interpreted as an operator 
in the (NS, R) sector and ip in the (R, NS) sector in the RNS formalism. Then, under 
the so-called (— l) Fi and (— 1) Fr transformations changing the signs of operators in the 
left-moving Ramond sector and those in the right-moving Ramond sector respectively, 
they transform as 

(_1)F £ . (7.1) 
(-1) Fr : if>->-ip, ip^ip. ( 7 - 2 ) 

In order for the matrix model action (12. ip to be invariant under the transformations, B 
and should transform as 

(-l) Fi : B^B, 0^-0, (7.3) 
(-l) Fi? : B ->B, 0^-0. (7.4) 

This indicates that B corresponds to an operator in the (NS, NS) sector, and in the (R, 
R) sector. In this interpretation, the Ramond sector is considered to be responsible for 
the logarithmic singular behavior seen in correlators among operators (12. lip . 

In fact, two-dimensional type II superstring theory with the identical target space, 
(<f,x) G (Liouville direction, S 1 with the self-dual radius), is constructed in [2U ESI [261 
127] . The (holomorphic) energy-momentum tensor on the string world-sheet except ghosts' 
part is 

Tm = ~\{dxf - - \{d V f + |<9V - \i>Mi (7.5) 

with Q = 2. ip x and if)g are superpartners of x and ip, respectively. Target-space super- 
currents in the type IIA theory 

q,(z) = e~^ zS) ~^ H ^~ ix{ - z \ q-(z) = e~^^ + i S ^ +i ^ (7.6) 

can exist only for the S 1 target space of the self-dual radius, which comes from the GSO 
projection implemented by imposing locality of the supercurrents with operators of the 
theory. (0) is the holomorphic (anti- holomorphic) bosonized superconformal ghost, and 
the fermions are bosonized as 

i)t ± iip x = V2e TiH , fa ± i<$ w = V2e TiR . (7.7) 

Then the supercharges 

Q + =l^-. q + (z), Q-=l^~. q-{z) (7.8) 



2ni " / 2ni 
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are nilpotent Q 2 + = Q 2 _ = 0, which indeed matches the supercharges Q and Q in the 
matrix model in (12. 2p and ( 12. 3|) . 

The spectrum except special massive states is represented by the NS "tachyon" vertex 
operator (in (—1) picture): 

rp^ _ e ~4>+ikx+p t <p rp_ _ e -4>+ikx+p t (p 

and by the R vertex operator (in (— |) picture): 

with e, e = ±1. Locality with the super currents, mutual locality, superconformal invari- 
ance (including the Dirac equation constraint) and the level matching condition determine 
physical vertex operators. As discussed in [26], there are two consistent sets of physical 
vertex operators - "momentum background" and "winding background" . Let us consider 
the "winding background" 0. The physical spectrum is given by 



(NS, NS) : 


Tk T-k 




(R+, R-) : 


Vfc,+i V-k,-i 


(k= 2> 2'"')' 


(R-, R+) : 


V-k,-l ^4,+l 


(fc = 0,1,2,---) 


(NS, R-) : 


T-k V-k,-i 


(k = -, -,-••)> 


(R+, NS) : 


Vh,+1 Tk 


(k= 2> 2'" -) ' 



(7.11) 

where we take a branch of pe = 1 — \ k\ satisfying the locality bound [40] . We can see that 
the vertex operators 



Vi +1 V_i T_iV_i _ 1: Vi +1 Ti, T 1T1 (7.12) 

2> ' 2> 2 2' 2' ' 2 2 2 

form a quartet under Q + and (Q-: 

v h + i = r_i VLi,-!, VL|,_ a } = o, 

{g + ,Vi i+1 f | } = T_if | , [g + ,r_ifi] = o, (7.13) 

{g_,r_i vii _!} = r 1T1, [q_,t Lift] = 0. (7.14) 



6 We can repeat the parallel argument for "momentum background" in the type IIB theory, which is 
equivalent to the "winding background" in the type IIA theory through T-duality with respect to the S 1 
direction. 

7 We here assume that Q + commutes with and anti-commutes with Vj, g , and that Q_ commutes 
with Tfc and anti-commutes with Vk, e - It is plausible from the statistics in the target space. In ref. |28) . 
we introduce cocycle factors to the vertex operators in order to realize the (anti-)commutation properties. 



25 



Notice that (I7.13P and (I7.14p are isomorphic to (12.21) and (12. 3p . respectively. It leads 
to correspondence of single-trace operators in the matrix model to integrated vertex op- 
erators in the type IIA theory: 




(7.15) 



which is consistent with the identification in ( I7.ip -( r7~4"j) . Furthermore, it is natural to 
extend (I7.15P to case of higher k{= 1, 2, ■ • • ) as 



Since the "tachyons" of the negative winding J d 2 zT_ k _i(z)T k+ i_{z) (k = 0,1,2,-..) 
are invariant under Q + and they are expected to be mapped to {^tr (— iB) k+1 } 
(k = 0, 1, 2, ■ ■ • ) perhaps with some mixing terms. We see in (I7.16P that the powers of 
matrices are interpreted as windings or momenta in the S l direction of the type IIA theory. 
Although such interpretation is not usual in matrix models for two-dimensional quantum 
gravity coupled to c < 1 matters, refs. [211 122] show that a positive power A; of a matrix 
variable in the Penner model correctly describe the "tachyons" with negative momentum 
— k in the c = 1 string on S 1 , which is in harmony with our interpretation. In [2~T| 122], the 
positive momentum "tachyons" are represented by introducing source terms of an external 
matrix via the Kontsevich-Miwa transformation in the Penner model. In turn, it is natural 
to expect in our case that the positive winding "tachyons" f cPzT_ k _i(z)T k+ i(z) [k = 
— 1, —2, ■ ■ ■ ) in the type IIA theory are expressed in a similar manner in the matrix model. 

Thus, from the argument here based on symmetry properties and spectrum of the 
matrix model and the type IIA theory, we have a plausible reason to expect correspondence 
between them. The vertex operators in the (R— , R+) sector are singlets under the target- 
space supersymmetries Q + and Q-, and seem to have no counterparts in the matrix model. 
As we have seen in sections El E] and El one- and three-point functions of $2fc+i and two- 
point functions of fermions ^2k+i, ^2fc+i are nonvanishing. From the point of view of 
our correspondence (I7.16p . it implies that correlators of operators with nonzero Ramond 
charges do not vanish, and hence the matrix model should correspond to the type IIA 
theory on a background of (R— , R+) operators. Furthermore (1A.25|) suggests that the 




+ (mixing) 



+ (mixing) <^ 



+ (mixing) 




(7.16) 
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torus free energy of the matrix model Fi in (I2.10p vanishes. Although it is different from 
the results of calculations in the string theory on the trivial background [2HE2], we can 
find a consistent (R— , R+) background in the type IIA theory which gives vanishing 
torus partition function in agreement with the result of the matrix model. Details will be 
presented in ref. [28] . 

So far we have seen the correspondence between our matrix model (12.1 j) and the 
two-dimensional type IIA theory with an RR background at the level of symmetries 
and spectrum. This correspondence is in fact more than that and quantitative. In the 
next paper |28j, we will calculate various amplitudes in the type IIA theory and find 
surprisingly that they also exhibit powers of logarithm as in the amplitudes in the matrix 
model presented in this paper. Thus we will directly see the correspondence quantitatively 
between amplitudes in both sides in f )7.16p . 

Our interpretation of the matrix model as superstrings does not rely on the picture 
of unstable D-branes discussed in [2]-[llj. It will be interesting to consider a relation of 
the matrix model with D-branes. In fact, it is quite likely that our matrix model would 
be interpreted as an effective action of D-branes and that the correspondence to the IIA 
theory could be a consequence of the open-closed duality. For the aim, ref. (13] would be 
helpful, which discusses FZZT branes HI] in the Kontsevich-Penner model. 
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A Analysis via Schwinger-Dyson equations 

In this appendix we analyze correlation functions of even-power operators of 4> via Schwinger- 
Dyson equations. For this purpose let us consider the Schwinger-Dyson equations 







.N N / \N 

k=0 x ' N 



/i_ tr (0-+i(^_^))\ (A.l) 
, , Y /ltr0"-tr/-t 

k=0 x 



= -IL ( i-tr 6 n+m )+y (^ti 4> n -^tr 6 k i-tr ^ ? 



which are derived from 

JV 2 



(A.2) 



(A.3) 



= fd N2 Bd N2 <t>d N2 i>d N2 TP^-^[tT(<t> m+ H a )e- s ], 

N 2 

= J d N2 B d N2 <j> d N2 iP d N2 4> J2 ^ m+1 n e~ S ] , (A.4) 

respectively. Here, is expanded by a basis of N x N hermitian matrices {t a } (a = 

l,---,iV 2 )as0 = E«=i0^ a 

By diagonalizing <\> and integrating out fermions, -^tr (0 m+1 (ipip — ipip)) in ( lA.ip and 



m+l , \ m + l 



( 1A.2j) becomes —-^2 J2ij=i * \ + \- — > where Aj (z = 1, • • • , N) are the eigenvalues of 0. 
Note that when m is even, this is further reduced to the polynomial: 

* J fc=0 

Thus, for m = 2£ (£ = 0, 1, 2, • • • ), (lA.lj) and (1A.2j) can be rewritten as 



£ (^tr0 2 * ^tr0 2 ^)) = i (ltr(50^ 2 )) , (A.6) 

k=0 



and 



1 ™ / l tr 0" +2 A +y( -tvr -tr 2fc ltr0 2 ^ fc ) 
= z/ltr^ltr(S^ 2 )^ (A.7) 
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respectively. 

In terms of the resolvent for <b 2 : 



R 2 ( z ) = —tT 1 



N z- 



eq. (1A.6j) is expressed as 



z (R 2 (z) 2 ) = (z 2 - fx 2 z) (R 2 (z)) - z + f x 2 -C, C = ( -tr 2 



1 



Also, f]A.7|) for n even is as 



(A. 



(A.9) 



with 



— d w (wD z ,(z, w) (R 2 (z'))) + z (R 2 (z) 2 R 2 (w)) 

{z 2 - fi 2 z) (R 2 (z) R 2 (w)) + {-z + fx 2 ) (R 2 (w)) ~ ( j^trn 2 R 2 ( ,r 



D Az , W )f(z>)= f{z) - f{w) 
z — w 



(A.10) 
(A.ll) 



A.l Solution of KKM and ( RTTUj ) 

Eqs. (1A.9|) , flA.lOj) can be solved by the genus expansion. By writing the genus expansion 
of C as in fl23|) : 



c = c + ±c 1 + ±c 2 + 



the O(N ) part of flAT9|) reads 



z ((R 2 (z)) f = (z 2 - fi 2 z) (R 2 {z)) - z + /i 2 - C . 



Its solution that does not blow up as z — > oo is 



(W)o = 2 



z-fx- \ (z - /j, 2 ) 2 - 4 + 



4(//2 - C ) 



(A.12) 



(A.13) 



(A.14) 



Assuming that (R 2 (z)) should be analytic except one cut on the positive real axis [a 2 , b 2 } 
(0 < a 2 < b 2 ) determines C as 

C = fx 2 . (A.15) 

Thus, we obtain 



z-n 2 - ^{z-a 2 ){z-V) 



(R2(Z)) = ~ 

a 2 = -2 + fx 2 , b 2 = 2 + /x 2 



(A.16) 
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which coincides with the result from the saddle point equation. Note that it holds for 
general filling fractions (u + , vJ). 
From the 0(N~ 2 ) part of (IP]) : 



(W%,0 " V / (Z-* 2 )(Z-V)(R2(Z)) 1 = - 

and the £>(A^ 2 ) part of 

d w w) (R 2 (z')) ) + z (R 2 (w)) {R 2 {zf) CQ 



Ci 

z 



-zy/(z-a*)(z-P) (R 2 (z) R 2 (w)) Ci0 + (R 2 (w)) (R^z)), 
1 



N 



tr<ff R 2 (w] 



Ci (R 2 (w)) 



C,0 



for both of which we have used (1A.16|) . we obtain 



d w (wD z ,(z,w) (R 2 (z')) ) - z^{z-a 2 ){z-b 2 ) (R 2 {z) R 2 (w)) 
= -/.Ltr^K 



CO 



Applying 



^)(^H)o = (^r0 2 J R 2 K 



and f TAIBj) to (IA~T9|) leads to 



(R 2 (z) R 2 {w)) 



Co 



A(z-w) 



(z — a 2 )(w — b 2 ) \ (z — b 2 ){w — a 2 ) 



(z — b 2 )(w — a 2 



[z-a 2 ){w-b 2 ) 



- 2 



It is easy to confirm that this equation yields (14. 13|) for p, q even. 
In particular, 

(M*) 2 ) Cfi = ( z _ a2 y( z _ b2 y> 
and requiring the regularity of (R 2 (z)) 1 at z = in ( 1A.17|) gives 

1 



d = o. 



(z- a 2 fl 2 (z-b 2 f/ 2 ' 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



(A.21) 



(A.22) 

(A.23) 
(A.24) 



Eq. (|A.24j) means that the torus free energy is a constant independent of fi 2 . We explicitly 
computed ^tr.B n ) (n = 1,2,3,4) from the large- z expansion of (1A.23[) . and obtained 
the result suggesting supersymmetry invariance at the torus topology: 



— tiB n 
N 



(n = 1,2,3,4). 



(A.25) 
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In ref. [5] we find strong evidences that the supersymmetry is preserved at the sphere 
topology in the asymmetric one-cut phase and in the two-cut phase both corresponding 
to (IA.16I) . It is therefore interesting if the supersymmetry remains intact at the torus 
topology. 



B Dependence on filling fractions of (Y[ i j^tr 



C 



We here discuss dependence on filling fractions (z/ + , vJ) of the .fT-point correlation func- 
tions (j[f =l itr (f) ni } with K — 1,2, 3. 

Let us express as Ze v „_} contribution to the partition function (12.41) from the config- 
urations where the first v + N eigenvalues of the matrix are around the right minimum 
and the remaining z/_iV around the left minimum: 



Xi = ji + Xi for i 6 I = {1, ■ ■ ■ , u + N}, 
Xj =-fi + Xj for j G J = {v + N + !,••• ,N} 



(B.l) 



with Xi and Xj small fluctuations. Then, 



N 

E 



(is + N)\(v_N)\ {v+ ' v - ] 



(B.2) 



u+N=0 

is valid at least to all orders in perturbation with respect to the small fluctuations, and 

N 

Z(y+j/-.) — v \ / (II | 11'' 'f > ," /x j ' 

2 



c N j (\[dxAX[w\ii + x i )X[w'\-ii + x ] 



vfc=l 

X 



el jeJ 

n (w'(i2+x t )-w\fi+x l 

i>i', 

x Y[ (w'i-fji + Xj) - W'i-fi + Xj> 

j>j',j,j'eJ 

X 



iei,je.J 
x exp 



II (w'(-fx + Xj)-W'( f j, + X i )^ 



(B.3) 



with W'(x) = x 2 — /i 2 , W"(x) = 2x. Cn is a numerical factor depending only on iV [3"T] . 
The integrals in (1B.3P are calculated in perturbation with respect to the small fluctuations, 
namely in the -^-expansion. By flipping the sign of Xj (j e J), it is easy to see 



(B.4) 
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Let us write correlators evaluated under the configurations (IB.ip as (•y +,v \ In 
the remaining of this appendix, we find a relation between ( Yli=i Ji^ Y ni ) an d 

r£i£tr0«*) for if = 1,2, 3. 

B.l K = l 

Applying the same argument of flipping the sign to (±tr(f> n ) iv+,1 '- ) leads to 



jr**")^ = (i* X> + ~ Xi)n ) + ( ~ 1)n ( ^ 5> + ~ Xj)n ) • (R5) 

Since there is a permutation symmetry of Ai, ■ • • , Ajy in (-}^ , we can write the r.h.s. as 

/ ~ \ (i.o) / 1 \ (1 ' 0) 

(u + + (-l)V.) ((ai + AO") = (*/ + + (-l)V.) ^tr^ • (B.6) 



(i.o) / _ \ (i.o) 



Thus 



I v («'+,«'-) / j \ (1,0) 



tr^M = (i/ + + (-l)V_W-tr^) (B.7) 



Namely, (^tr <p n ^ v+,v ^ is independent of the filling fractions for n even, and proportional 
to (u + — u-) for n odd. Note that the relation (IB .71) holds at each order of the -^.-expansion. 



B.2 K = 2 

After the sign flip, we have 

I I \ ("+>"-) 



-tr0--tr0- 



c 

(1,0) 



/ 1 1 \ (1 ' 0) 

+(-i) ni+na ( jj z> + at Eo* + ) ( B -9) 

\ iieJ 7' 2 eJ / 



(1,0) 



( ^ EC" + ^) ni ^7 £(*« + ^) n2 ) ( R1 °) 



A" A 



(1,0) 



( a7 £(<" + ~ X ^ ni 1j £0" + ^ > ( B - n ) 



! v (1,0) / 1 \ (i,o) 



A ^— ' ' x/ A 

(z/ + + (-irV_)(z/ + + (-l)"V_)^tr0" 1 \ (^tr0« 2 ) . (B.12) 
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Use of the permutation symmetry recasts the first term in the r.h.s. as 

(1,0) 



\"2 



(HD = ^((/i + A 1 ) ni (^ + A 2 ) ? 

^ ( { (jm + - (/i + A 2 r } { (/x + AO- - (/i + a 2 )- 2 } )^ 0) • 

(B.13) 



+ 



Notice that the connected part is taken in the second term because the disconnected part 
vanishes due to the permutation symmetry. Since (z/ + ,z/_) = (1,0) case in the above 
argument yields 



/ 1 1 \ u,u; / \ (i,o) 

\;v tr<r iv tr r / = + Al)ni (yU + Aa) / 

({ {fi + ~ Xi)ni - ^ + ~ X2)ni } { {fi + ~ Xi)n2 - ^ + X2)n2 })T ' (B - i4) 



we obtain 



/I i \ (i,o) 



(M) = ^+( — tr0 ni — tr0" 2 



+^ ( { + AO" 1 - (A* + A 2 ) ni } { (/x + AO" 2 - (/i + A 2 )" 2 } ) 



(1,0) 

c* 



(B.15) 



Repeating the same argument for (1B.9|) . (IB.lOp and (IB. lip , and putting the results 
together with the last term flB.12j) . we end up with 



-tr^-tr^ 

= (u + + (-i)-v_) (u + + (-i)" 2 ^) fair ^ tr ^ 2 )^ 0) 
+(i-Hn(i-(-in^ 

x ({(/x + A0 ni + A 2 ) ni } {(// + AO" 2 - (/* + A 2 ) n2 })| ; 1 ' 0) • (B.16) 
This formula tells us that 

i i \ I i i \ (1,0) 

-tr^ ^tr^\ = /-tr0- ^tr^\ (B.17) 



for ni and 712 even, and 



i i \ (f+f"-) /i i \ (i,°) 

-trc^ -tr0« 2 \ = („ + - »,_) / -tr<T -tr<T \ (B.18) 

i C > ' c 
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for one of n\ and n 2 odd. These are valid to all orders in the -^-expansion. For n\ and n 2 
odd, the second term in the r.h.s. of (IB. 16[) appears to be nonvanishing. However, it is at 
most of 0(N~ 3 ) and negligible compared to the leading contribution of 0(N~ 2 ). Thus 
we have shown that 



(?+,"-) I 1 i \ (i.o) 



for rii and n 2 odd. 



(B.19) 



B.3 K = 3 

The argument parallel to the previous subsection leads to the result 



Vt=l 



riK+(-ir-)}(n^ 



■tr. 



(i.o) 



c 



' 2N 



C k. i=l J \i=l 

(l - (-ir ) (i - (-ID iy+ + (-i) n v_) 

{ (// + AO" 1 - (/x + A 2 ) ni } { (// + AO" 2 - (// + A 2 )" 2 } (// + A 3 ) 
+ (two terms from cyclic permutations of (tii,ti 2 , n 3 )) 



(i.o) 



c 



, 2z7 +^ ft \ 

+ „ w f{ni,n 2 ,n 3 ) 



3N 2 



{ Oi + A x ) ni - (a* + A 2 ) ni | |(/i + Ai) n2 - (/i + A 2 )" 2 } 

x {(/i + Axr + (// + A 2 ) ns - 20 + A 3 r' x lJJ " 



c 



+ (two terms from cyclic permutations of (ni, n 2 , n 3 )) , (B.20) 



where 



f{n 1 ,n 2 ,n 3 ) 



(two or all of n\, n 2 , n 3 are even), 

1 (one of rii, n 2 , n 3 is even), 
z/ + — i/_ (all of ni, ^2> n 3 are odd). 



(B.21) 



From (IB .2 Op . we see that 



iHtr^ 



,i=i 



e 



v£=l 



(1,0) 



c 



(B.22) 



for Tii, n 2 and ra 3 all even, and 



vi=l 



ri> 

vi=l 



(1,0) 



(B.23) 
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for one of n\, n 2 and n% odd. These hold to all orders in the -^-expansion. For the 
other cases, although the second and third terms in the r.h.s. of (1B.20I) appear to remain 
nonzero, they are at most of 0(N~ 5 ). So, by looking at the leading order of 0(N~ 4 ), 

3 \ ("+•"-) / 3 \ ( X '°) 

II ^ tr ) =("+- *-)* ( II ^ tr ^ ) (B.24) 

i* =1 / C,0 \* =1 / C,Q 

is shown. Here, fl is the number of odd integers in {711,77,2,773}. 



rK 1 

U=l N 

general K is proportional to (v+ — vJf- with (J the number of odd-power operators of 



From the results obtained in this section, it is expected that { Yl i=1 j^ti <p ni j for 



C Useful formulas 

The semi-circle eigenvalue distribution yl — y 2 [y € [—1,1]) appearing in the Gaussian 
one-matrix model S = Ntr (2(f) 2 ) satisfies the saddle point equation 



I dy\'l — y 2 P = irx for x e [— 1,1], 

7-i x-y 



(C.l) 



from which we obtain 



1 1 

dy V 1 - y 2 p y 

1 x-y 



x 



y i dy V^fP^ ~ y i dy = n ( x2 -\)- ( C - 2 ) 



Next, let us compute 



/ dy^/l-y 2 ?—?—, x, we [-1,1]. 

J -1 x-y u-y 



(C.3) 



We express the principal values as 

p . 1 1 ( 1 1 



x — y 2 \x — y + ie x — y — ie 
1 1/1 1 



u — y 2 \u — y + ie' u — y — ie' 
with e and e' infinitesimal, and rewrite P— !— P^— as 

' x— y u—y 



(C.4) 



x—y u—y 4 



11 11 

+ 



x — y + ieu — y + ie' x — y — ie u — y — ie' 

11 11 
+ ; + 



x — y + ieu — y — ie' x — y — ie u — y + ie' 
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(C.5) 



We recast the first term as 
1 1 



x — y + ie u — y + ie' 
1 



x — y + ie u — y + ie' J u — x + i(e' — e) 
1 



- in5(x — y) — P h in5(u — y) 

x — y u — y 



u — x + %{e' — e] 



(C.6) 



and carry out the integral in the l.h.s. of (10.3|) by using ( IC.ip . After doing in the same 
manner for the other terms in (IC.5I) . we arrive at 



7T , - 

TC + i — v 1 — X 2 
4 



1 



+ 



1 



u — i 4" i(e' — e) u — x + i{—e' + e 
1 1 



+ 



u — x + i(e' — e) u — x + i(—e' + e) 
1 1 



+ 



u — x + i(—e' — e) u — x + i(e' + e 



(C.7) 



Here, since (IC.7I) becomes equal to — tc + tx 2 \/\ — u 2 5{u — x) in the limit e, e' — >■ with 
either e > e' > or e' > e > 0, we conclude that 



/i ^ ^ 
dy Jl-y 2 Y P = -it + 7rVl - u 2 5(u - x) 
. x x-y u-y 



(C.8) 



for x,u E [—1, 1]. 
Finally, we show 



' ,P^^ = for xe[-l,l]. 



dy 

i Jl-y 2 x ~y 



(C.9) 



The l.h.s. is written as 



(l.h.s.) 



1 1 



2 1 - x 
1 1 



21 + x 



/ dy^y~ 2 (P— - -L-) 
J-i V. x-y l-yj 

J^dyVT^(p^- 



y -i-y 



(CIO) 



Applying (IC.lj) leads to 

(l.h.s.; 



ii .11. 

\HX — IT) H (71X + 7TJ = 0. 

21 -x V ; 21 + x V ; 



(C.ll) 
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